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ABSTRACT 

Halo Abundance Matching has been used to construct a one-parameter mapping be¬ 
tween galaxies and dark matter haloes by assuming that halo mass and galaxy lu¬ 
minosity (or stellar mass) are monotonically related. While this approach has been 
reasonably successful, it is known that galaxies must be described by at least two 
parameters, as can be seen from the two-parameter Fundamental Plane on which 
massive early-type galaxies lie. In this paper, we derive a connection between initial 
dark matter density perturbations in the early universe and present-day virialized dark 
matter haloes by assuming simple spherical collapse combined with conservation of 
mass and energy. We find that z = 0 halo concentration, or alternatively the inner 
slope of the halo density profile a, is monotonically and positively correlated with 
the collapse redshift of the halo. This is qualitatively similar to the findings of some 
previous works based on numerical simulations, with which we compare our results. 
We then describe how the halo mass and concentration (or inner slope a) can be used 
as two halo parameters in combination with two parameters of early-type galaxies to 
create an improved abundance matching scheme. 

Key words: galaxies: haloes — galaxies: evolution — galaxies: formation — dark 
matter — cosmology: theory 


1 INTRODUCTION 

Halo abundance matching is one of several methods 
used to link galaxies with dark matter haloes. It uses 
the simple assumption that galaxy luminosity (or stel¬ 
lar mass) and halo mass are monotonically related, such 
that more luminous galaxies reside in more massive haloes, 
to match observed magnitude-limited samples of galaxies 
to dark mat ter halo merger trees fr om dark matter-only 
simulations llVa le fc Ostrikerl [^04al: iK ravtsov et al.l 120041 : 


IVale fc Ostrikerf 2004bl: Guo et al. 201 Oh. Despite the 

plicity of its underlying assumption, abundance match¬ 
ing is able to reproduce with surprising accuracy vari¬ 
ous measures related to the observed physical distribution 
of galaxies, such as the luminosity functions of different 
cosmic environments, t he occupation numbers of haloes 
I Vale fc Ostrikerl 2004a), galaxy autocorrelation function s 
I Conroy et al.l 120061 : IGuo et al.l l2010l: iN uza, et a l.l l2013ll. 
and g alaxy-galaxy lensing ( Hearin et al.lkOl^ Fosalba et al.l 

Being able to match large samples of galaxies and 
haloes without the use of complex semi-analytic or numeri- 
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cal hydrodynamic modeling has allowed for a simple probe 
of the connection between galaxies and their dark mat¬ 
ter haloes. Abundance matching has been used to obtain 
various statistical relationships between galaxies and their 
host halos, suc h as the luminosity-halo ma.ss and stellar-halo 
mass relations (TVale &: Ostrikerl 2004bl: Shankar et al. 200d 


Conroy fc Wechsler 20091: Behroozi et al. l2010l: Guo et al 


201C : IWake et al.ll201lF Leauthaud et al.ll2012l ; lMoster etTd 

2013r) . the relationship between galaxy optical circular ve¬ 
locities and the circular ve locities implied from their host 
haloes llDutton et al.l l2010ll. the halo baryonic mass func¬ 


tion (iBaldrv et al.l 120081 1. and t he relation between central 
black hole mass and halo mass llShankar et al.lBoOfill . It has 
also been used to examine the dark matter haloes hostin g 
certain types of galaxies, such as quasars (ICrotonI l2009ll . 
Furthermore, abundance matching also allows the possibil¬ 
ity of assigning observed galaxies at different redshifts to 
simulated haloes whose mass growth and merger history 
are known, allowing one to track the evolution of galaxies 
through time. This has been used to study the fate of satel¬ 
lite galaxies in clusters (IConrov et al.l2 007ll. the frequen cy of 
gas-rich versus gas-poor mergers I Stewart et al.ll2009bl l . the 
accuracy of obser vational indicators of halo mergers such as 
close pair counts llStewart et al.|[2^009al l . the evolution of the 
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stellar-halo mass relation and velocity dispersion-halo mass 
relation and the implications for galaxy- halo co-evolution 
(jFirmani fc Avila-Reesell201f)l : IChae 1201 ih. an d the growth 
of Brightest Cluster Galaxies ( Lanorte et al.l[2013l l. among 
many others. 

However, a scheme that treats luminosity as the sole 
important property of a galaxy cannot be entirely cor¬ 
rect. Galaxies are described by at least two parameters, 
as is apparent for moderate mass systems from the bi- 
modal distribution of galaxy colors at fixed lumin osity (e.g., 
IStrateva et al.l[200ll : iBlanton fc Moustak^l2009l l. This im¬ 
plies that at least one other parameter of the galaxy’s halo, 
aside from its mass, must be relevant to the evolution of the 
galaxy. 

Some recent work has explored the addition of an¬ 
other halo and galaxy parameter to abundance match¬ 
ing at fixed stellar and halo mass (iHearin fc WatsCTill2013l : 
iHearin et al.l l2013l : IWatson et al.l I20l4l : the authors refer 
to this scheme as “con ditio nal abundan c e ma tching.” In 
IHearin fc WatsonI ll2013tl and IHearin et all ll2013l ). the colors 
of galaxies are matched to a proxy for the halo age. The lat¬ 
ter is parametrized as the earliest of three times: when the 
main halo progenitor mass exceeded 10^^ Mq, when the halo 
became a subhalo, or when the halo transitioned from fast 
to slow dark matter accretion, w hich is computed dire ctly 
from the halo concentration as in IWechsler et aP (l2002tl . In 
practice, the last, concentration-based, age parameter is the 
one used for all but the most massive galaxies. The authors 
find that their method is able to match a number of observ¬ 
ables for galaxies separated into blue and red colors, includ- 
ing clustering statist ics and the galaxy-galaxy leasing signal. 
IWatson et al.l ll2014l ) apply the same method, but matching 


star formation rate (SFR) instead of galaxy color to halo age. 
They also report good agreement with observations, largely 
due to the correlation between SFR and galaxy color. 

Even considering only massive early-type galaxies, one 
parameter is inadequate for predicting all of their proper¬ 
ties. It has been known for some time that, while a rough 
one-parameter relation exists relating all variables to the ve- 
locity dispersion (the well-known “Faber-Jackson relation;” 
iFaber fc Jacks(^ll976l l. a two-dimensional parametrization 
called the “Fundamental Plan e” (FP) offers a superior de¬ 
scription (IPressler et aLlflOSTl l. An ev olving FP has been 
detec ted from z = 0 out to 2 ~ 2 (Ivan de Sande et al.l 
I 2 OI 4 II . Other properties of early-type galaxies, such as ob¬ 
served galaxy color, as well as modeled stellar population 
ages and metal abundances, have b een found to be high ly 
correlated with the F P parameters (iGraves et al.l l2009al la: 
iGraves fc Faberil201(J) . This implies that the properties of 
early-type galaxies may be well-described by two parame¬ 
ters, making them a good sample on which to test a two- 
parameter matching scheme. 

There are also other reasons why a two-parameter 
matching scheme may work best for massive early-type 
galaxies. The disk-to-bulge ratios of spiral gal axies are likely 
to be dependent on their environment (e.g., iHopkins et al.l 
I2OO9II . Also, galaxy color in star-forming galaxies will not be 
well-correlated with stellar age because a recent small burst 
of star formation can make a galaxy significantly bluer while 
only slightly changing the mean stellar population age. 

One halo parameter of physical interest is the col¬ 
lapse time of the dark matter halo. That the properties of 


galaxies hosted by dark matter haloes could be affected by 
variations in collapse time at fixed halo mass is an idea 
that has appeared in a number of previous works (e.g., 
iBlumenthal et aI1ll984 l. In particular, since the galaxy color 
of early-types should be well-correlated with the age of the 
stars in the galaxy (as well as their metallicity), it is possi¬ 
ble that the galaxy color is also correlated with some mea¬ 
sure of the collapse time of the host halo, which would de¬ 
termine when gas could collapse and form stars. There are 
also other galaxy parameters that are likely to be correlated 
with the collapse time of the host halo at fixed mass, such as 
the metallicity and stellar mass-to-light ratio. A present-day 
halo property correlated with the halo collapse time could 
then be matched with a present-day galaxy property as a 
second set of abundance matching parameters. 

In this paper, we use a simple spherical collapse model 
to derive a present-day halo parameter that is a proxy for the 
halo collapse time. We adopt two different fitting functions 
for the 2 = 0 halo density profile, and also consider physical 
parameters that are independent of the function used to 
fit the halo profile. We will make use of the parameters we 
have found as part of a matching scheme to the Fundamental 
Plane in a future paper. 

We first derive a proxy for halo collapse time using 
simple spherical collapse model in 321 We show the results 
of this model in El We compare our resul ts with previ¬ 
ous parametrizations of halo “formation time” dBullock et al.l 
I2OO1I : IWechsler et al.l I2OO2I : IZhao et al.l l2009ll derived from 
dark matter simulations in El Finally, we describe a scheme 
for halo matching of the Fundamental Plane of elliptical 
galaxies (a two-parameter distribution), which is not obvi¬ 
ously dependent on environment, that we plan to expand on 
in a future paper in 311 

In all parts of this paper we assume a cosmology consis¬ 
tent with the WMAP nine-year results plus external GMB, 
BAG, and Ho measurements (iHinshaw et al.l l2013l . Table 
4); thus we take Ha,o = 0.71, Hm,o = 0.29, and Hq ~ 69 
km/s/Mpc. 


2 SPHERICAL COLLAPSE MODEL FOR 
HALO COLLAPSE TIME 

We would like to choose some property of dark matter haloes 
that is a good proxy for the halo collapse time and can be 
easily measured in dark matter-only simulations. We ap¬ 
proximate the collapse time of a dark matter perturbation 
early in the universe as twice the turnaround time in sim¬ 
ple spherical model s of collapse. W e use a somewhat similar 
method to that of iRubin fc LoebI (I 2 OI 3 I I, who give equa¬ 
tions for calculating the virialization density Ac for arbitrary 
pre-collapse and post-collapse density profiles by assuming 
mass and energy conservation. Our method is similar, ex¬ 
cept that we match the initial and final profiles within their 
turnaround radii at 2 = 0. 

For the purposes of creating a two-parameter match¬ 
ing scheme, we choose final halo profiles that are described 
by two parameters and that are commonly used to fit nu¬ 
merically simulated haloes—namely, an NFW profile with 
parameters M 200 and c, and a generalized NFW profile with 
varying mass M 200 and inner slope a, with fixed c = 5 (see 
32 : 211 . We match these final profiles such that their mass M 
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and energy E within the shell that turns around at « = 0 
is the same same as that for chosen arbitrary initial profiles 
(e.g., tophat or Gaussian). While this model maps final pro¬ 
files to initial profiles with two parameters uniquely (details 
below), the corresponding final profile is not actually the 
profile that the initial profile would evolve to, as evidenced 
by the fact that profiles of different shapes can be matched 
to the initial profiles this way. Rather, we choose final pro¬ 
files that are used as approximations for a variety of dark 
matter halo profile shapes. 


2.1 Initial Profile 


Here we review the equations for the ev olution of t he ini tial 
profile; the analysis is similar to that in I Mo et al.l (l2010h . 

We begin with some chosen initial overdensity profile 
Pipi) at arbitrarily chosen initial time ti. We assume the 
profile will tend to the mean matter density of the universe 
at that time, p(ti), for large radii. As long as the density 
is decreasing or constant with increasing radius, there will 
be no shell crossing for shells that have not yet collapsed 
and we can treat them separately. It is assumed that going 
far enough back in time, the initial perturbation is entirely 
expanding, and none of its shells have yet turned around. 
Also, in a ACDM universe, a finite amount of mass will 
collapse in an infinite time, because for overdensities lower 
than some value, the shells expand forever due to the A 
term instead of collapse. In the Appendix to this paper, we 
present an exact derivation of this value for an initial profile 
assumed to be on the Hubble flow, which tends toward the 
solution having r = 0 at t = 0 for ti —>■ 0. 

The collapse of each shell enclosing mass M(< r) is 
governed by the following equation: 


cPr 

dP 


GM 2 

-5-f slA.orzor. 


( 1 ) 


The above equation integrated once becomes 
1 f dr y GM tlA.oHi 2 ^ 

2 UJ -= 


where £ is the specific energy of the shell. At the turnaround 
time tta of a given shell, this becomes 


£ 


GM 

Pta 


2 


2 


Tta, 


(3) 


where the turnaround radius rta is the maximum radius at¬ 
tained by each shell. Defining 

^ “ 2GM ’ '■ ’ 

we see that for r < 0 at rta we require ^ < 1/2. From this 
and Eqn.|2]we have a formula relating the radius of the shell 
and the time, as long as the shell has not yet turned around: 


Hot 



1 

X 


l + C.{x^ 



(5) 

Using Eqn.[5] the central collapse time (2 x tta) can be 
calculated for rta —>■ 0 (as long as the density profile does not 
have a central cusp), as can the time at which any fraction 
of the mass at 2 = 0 collapsed. 

For a selected initial density profile p{ri) at time U, we 


can find the energy of each shell from Eqn.O We first obtain 
rta for each r^. To do that we re-express ^ as 

, flA.oHg frtay 

^ yGpi{< ri)\n) ^ > 

and insert into Eqn. [5]for ti{ri). This can be solved numer¬ 
ically for rta as a function of ti and ri. 

By setting t in Equation [5] to the age of the universe 
at 2 = 0 (referred to here as to) we can obtain () of the 
shell that is turning around at 2 = 0. Then we can use 
this same equation to find the ratio of the original radius to 
the turnaround radius ri/rta of this shell by setting t = ti. 
This combined with Equation [U] gives us the ri and rta for 
the shell turning around at 2 = 0; we will designate these 
as ri^rnax and rta,max- We designate the mass and energy 
within ri,max at U (and rta,max at to) as Mtot and Etot, and 
these are what we will match to the mass and energy of the 
final profile. We also define the initial density within ri,max 
to be pi,max, such that 


^ f Mtot \ 

\4:/37Vpi,max J 


(7) 


We can find £ for each shell using Eqn. [S] The total 
energy of all the shells within the maximum radius will then 
be 


Etot 



£{rta)dM{rta). 


( 8 ) 


To obtain £{rta) for each shell at initial radius rt between 
ri = 0 and rt = ri,max, one can again insert Eqn. [6] into 
Eqn.[S]for t = ti, and solve for rta/ri. 

Scaling the size of the profile by the radius ri,max, so 
that the coordinate used is y = r/ri,max, Eqn. |3] gives 


^ = -^7rGpi(< y)—y^rlraax 
3 rta 


nA,oHi (rjaV ,0 (gx 

* I y ' i,max ? J 

n ' 


where rtapi is a function of y and also depends on 20 and 


Then, we can use Eqn. [U] to substitute for rt/rta in the 
equation for f; we obtain: 




Then the total energy is then given by 

ri 


/ 


Etot = / £dM = dTrr,; 




= 47rr“ 


1 

'2 V3 


£{y)pi{y)y dy 

2/3 

ttG) 


^ pV^i<y)piiy)^^^/^y*dy 

where using Eqn. 0 we then obtain 

Etot/M^^ = -^^^{HoG)^/^Hyx 

Pi,max 


( 10 ) 


( 11 ) 


^ pyi<y)piiy)^^^jS^y^dy. ( 12 ) 
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This equation holds for any initial density profile pipi) at 
any chosen Zi, as long as no part of the profile has yet col¬ 
lapsed. 

For the initial profile, we also want to consider the limit 
as Zi becomes large. Here 5(< y), wh ere pi{< y) = p{zi){l + 
5{< y)), approaches iIMo et al.ll2010l') : 


5 (< y) 


3 i + C(y) 

5 C/^y) 


^m,0 / 


{1 + Zi) 


(13) 


where corresponds to a given turnaround (or collapse) 
time. We want to consider the same overdensity profile shape 
at all times; i.e., 5{< y)ISmax does not vary with time, 
where Smax is the overdensity at ri^rnax- It is clear that for a 
shell with a given collapse time, 5 for that shell evolves with 
time; however, the above equation shows that for large Zi, a 
5(y)/Smax prohle taken to be constant with time is in fact 
also a constant profile in C,{y) —that is, collapse time at any 
fixed interior mass. 

Here pi{y) changes with «i, but for large Zi, pi{y) —>■ 
p{zi) as the perturbation 5 decreases like (1-|-Zi)“^ for fixed 
C (Eqn. [T3|, eliminating pi^max, pi{< y), and pi{y). This 
gives 


the mean matter density: 


p{r) 


_ Po _ 

r/rsil + r/r,)^ 


+ p{to)- 


(15) 


As noted above, for large values of the concentration c, 
oc c/(logc)^ for the NFW profile, which is not 
bounded. Therefore NFW profiles above some concentration 
cannot be the product of a simple spherical collapse model. 

(ii) While the NFW profile is the most common 
parametrization of dark matter halo density, others have 
found that haloes are equally or bette r described by pro¬ 
files t hat have a varyi ng inner slope (e.g.lSubramanian et al.l 
1200(f) . In particular, iRicotti et al.l 1 I 2 OO 7 I') have found that 
simulated dark matter haloes at virialization are equally 
well fit by either an NFW profile with varying c, or by a 
generalized NFW profile with fixed c = 5, and varying inner 
slope a. Thus we also match initial profiles to a dark matter 
density profile given by 


p(r) 


_Po_ 

(r/ra)“(l -I- rjrs)^~°‘ 


+ p{to), 


(16) 


IP /Ar5/3_ 3 ^^2/3nl/3 1 + Co S{<y) .ij.. 

S‘A,0 .1/3 / a y 

^ Jo ^max 


(14) 


where ("o is the value of C, for the shell that turns around at 
zo, in our case taken to be 2 = 0. 

Since ((0 is dependent only on the reference redshift zq 
that we select and we assume that the overdensity profile 
is a constant shape, this expression approaches a constant 
value as 2 i 00 . We see that is dependent only on 

the shape of the overdensity profile, which also determines 
the collapse times relative to the chosen reference redshift 
at high Zi. The lack of dependence on Zi for high Zi is due 
to the fact that A becomes unimportant at early times. 

It is also clear that for a constant or decreasing over¬ 
density profile, the integral on the right side of Equation 1 141 
is bounded between 1/5 and 1/2, meaning that the energy 
within a shell turning around at a given redshift must be 
bounded between two values. However, as described in the 
next section, the final (2 = 0) halo prohles we consider are 
common empirical fits to simulated haloes, and thus have 
unbounded possible values for . Therefore, not all 

conceivable final profiles are able to correspond to a possi¬ 
ble initial profile. 


2.2 Final Profile 

We would like to find a corresponding present-day halo 
density profile that has equal Mtot and Etot within the 
turnaround radius at 2 = 0 (t = to) as the initial profile. 
We note that a specific turnaround time defines a unique 
density within the turnaround radius at that time (Eqns. 
[4]and[5|, implying that for fixed Mtot, the final and initial 
profiles also have the same turnaround radius. 

We choose two different forms for the present-day halo, 
based on the fact that these shapes are commonly used to 
fit simulated haloes: 

(i) The entire profile is described by an NFW prohle over 


where Ts is related to M 200 via the hxed concentration c = 5. 
We refer to this prohle as an a model to avoid confusion with 
the NFW prohle. This prohle and the standard NFW have 
the same slope of —3 for large r, making their properties 
similar at large radii. As for the NFW prohle, is 

unbounded for large a. 

We assume haloes to be virialized within r 2 oo, the ra¬ 
dius within which the mean density is 200 times the critical 
density, and that the virial theorem can be used to hnd the 
energy within this region. Outside the virial radius, the pro¬ 
hle is collapsing, out to a radius rta at which tta = to- While 
a different radius could be chosen within which the prohle 
is virialized, we note that taking the viria l radius to be Voir 
as computed in iBrvan &: Normeml (Il998l l would produce a 
negligible difference in our results. 

The energy of the NFW prohle within the radius turn¬ 
ing around at 2 = 0 is the sum of the energy within the 
virial radius and the energy in the shells turning around. 
The potential energy within the virial radius is the sum of 
that from the matter and that from A: 

h r 

r ^'^00 1 

Ua = - -nA,oH^47vp{r)Pdr, (18) 

Jo 4 

and by the virial theorem, E{r 2 oo) = 0.5[/m + “JUa- Because 
p is much lower than the virial density, the fact that we take 
the density prohle to be an NFW or an a model in overden¬ 
sity has an insignihcant effect on the virial energy we hnd 
within the virial radius. Thus the energies we obtain within 
the virial radius for the NFW prohle are approximat e ly tho se 
given by equations 3.33 and 4.20 in I Rubin fc Loel^ ll2013l l. 

Outside the virial radius, the region between r 2 oo and 
the turnaround radius rta is collapsing. To determine the 
energy for this region we use Eqn. 1111 However, because 
the region is collapsing and not expanding, one must add 
the turnaround time to the time a shell has been collapsing 
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Figure 1. Correspondence between concentration c of NFW pro¬ 
files and inner slope of a models, for profiles at 2 : = 0 that have 
the same mass and energy within the turnaround radius. For 
comparison, gray p oints are from fits to simulated haloes from 
iRicotti et al.l <l2007l) . described in more detail in the text. 


after Ua, thus substituting Eqn.[5]with: 


Hoto = 


( c 

yi r 

/ dx 


[Jo [ 


1 + C(a:^ 


1 ) 


- 1/2 



i-l + C(x^-l) 

-1/2 

/ dx 


Jr/rta 

X 



(19) 


For a selected initial profile shape (e.g., a Gaussian) 
at a chosen initial time ti, the above steps will create a 
one-to-one mapping between initial profiles with mass Mtot 
and energy Etot and final profiles at 2 = 0 with the same 
mass and energy within the turnaround radius. Since both 
the NFW profiles and a model profiles are a two-parameter 
family, a unique combination of mass and energy values will 
correspond to a unique final profile of a given form. 


3 RESULTS 

Applying the above procedure, we are able to match final 
profiles to initial profiles. It is clear from Eqn. [T^] that the 
value at a fixed redshift is dependent only on the 

shape of the overdensity profile (and the mean cosmic den¬ 
sity at the chosen redshift), as it must be from dimensional 
arguments. For the final profiles we take to be at 2 = 0, NFW 
profiles and a models with the same mass and energy, this 
equation implies that there will be a relationship between c 
and a that is not mass-dependent. We show this relationship 
in Figure [ 1 ] More concentrated NFW profiles correspond to 
a model profiles with steeper inner slope. For comparison, 
we also show mean fits to the 40 most massive haloes from 
simulations at different redshifts from IRicotti et ^ ll2007t) . 
Our mapping using the energy and mass corresponds well 
to the match between a and c from direct fitting. This is 
as could be expected, since direct fitting ensures that the 
profiles will have similar shapes in the region r ~ r^, and 
both profiles have slope -3 at large radii, leading to similar 
energies at fixed mass for both profile shapes. 



1.98 

1.94 

1.88 

1.82 

1.73 

1.62 

1.46 

1.20 

0.72 


e 


Figure 2. Top Panel: The concentration c of final NFW profiles, 
or inner slope a of a models, for different central and half-mass 
collapse redshifts of initial Gaussian overdensity profiles at Zi = 
1000. The matching between initial and final profiles is described 
in the text. The 2 = 0 halo concentration (or alternatively, inner 
slope) is monotonically and positively correlated with the collapse 
redshift of any fraction of the 2 = 0 mass. 

Bottom Panel: The normalized squared maximum circular ve¬ 
locity, u^j,j,/(GA720o/A20o)i at 2 = 0 for an NFW profile with 
varying c, and an a model with varying inner slope a, versus the 
central collapse redshift of the matching initial Gaussian overden¬ 
sity profile with Zi = 1000. The range in c and a shown is the 
same as in the top panel. The similarity in v^^^l{GM2ooIE 200 ) 
for NFW profiles and a models over the range of physical interest 
is apparent. 


Similarly, Equation [12] implies that for the initial pro¬ 
files, E/M^!^ is also a function of only the shape of the 
overdensity profile and the mean density at the chosen ini¬ 
tial redshift Zi. Furthermore, as seen in Equation 1141 for a 
fixed overdensity shape, the value of approaches a 

constant value as Zi — >■ 00. However, the collapse time for a 
mass shell containing a certain fraction of the total mass is 
also a function of only the chosen initial redshift zi and the 
shape of the overdensity profile (Equations [5] and |6|) , and 
also approaches a constant value for fixed overdensity shape 
as 2i ^ 00 (Equation 1131) . Thus the value of c or a of the 
final profile will be a function of the collapse time of any cho¬ 
sen fraction of the mass for an initial profile of a fixed shape 
for fixed Zi. We show this correspondence in the top panel of 
Figure [ 2 ] in which we present the NFW concentration c (or 
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a model inner slope) versus the collapse time of the center 
of an initial Gaussian prohle, as well as the time for half 
the mass collapsed at z = 0 of the initial profile to collapse. 
For the initial profile we take Zi = 1000 and compute the 
exact collapse times as described in am however, this red- 
shift is large enough that the results will be similar to those 
calculated with the approximation of —>■ oo. We find that 
more concentrated haloes (or those with steeper a) have ear¬ 
lier collapse times, as might have been intuitively expected. 
Thus, either a ox c can be used as a proxy monotonically 
related to collapse epoch that is independent of halo mass. 

The results depend on the choice of initial profile. They 
also depend (weakly) on Zi, the initial redshift at which the 
profile is selected to be the shape of choice. This is because 
in general the shape of the profile evolves over time, so that 
the same profile at a later or earlier time does not follow 
the same functional form. Thus the dependence on Zi can 
also be seen as equivalent to a dependence on the shape of 
the profile at any given time. However, as seen in ^2.11 as 
Zi —>■ oo, the dependence on Zi disappears, and the value of 
depends only on the overdensity profile shape. 

Additionally, we present two tables showing the same 
values as the top panel of Figure [21 Tables [T] and [2] show 
the central and half mass collapse time, respectively, for 
NFW profiles with concentration c and mass M 2 oo. While 
the collapse times do not depend on M200, in bold we show 
the mass-co ncentration relation for sim ulated dark matter 
haloes from iDiemer Kravtsov! ll2014l ') , in an observation- 
ally normalized CDM universe within which both typical 
values of c and M200 are functions of collapse e poch. We 
select the results from iDiemer Sz Kravtsov! (l2014l l that use 
the same cosmology as we assume throughout the paper 
llHinshaw et al.l l2013tl . The range in c at fixed mass rep¬ 
resents their reported one-sigma scatter of 0.16 dex. By na¬ 
ture, our calculation is done for an arbitrary initial profile, 
so we do not assume cosmological initial conditions. Thus we 
must take the mass-concentration relation from elsewhere. 
For the assumptions about the initial profile we have made 
above, the bold values represent the expected scatter in the 
central and half-mass collapse redshifts for NFW profiles of 
a given virial mass. The ratio of the actual concentration of 
a halo to the mean concentration of haloes of its virial mass, 
c/c{M2oo), could be used as a parameter in a two-parameter 
matching scheme between haloes and galaxies that would be 
correlated with the halo collapse time and potentially with 
galaxy properties at fixed mass, as described in our discus¬ 
sion of future work in ^ 

The bottom panel of Figure [2] shows another parameter 
of the final profile, the normalized squared maximum circu¬ 
lar velocity v^ax/{GM 200 /R20o)j versus the central collapse 
redshift of the same initial profile as in the top panel. Again, 
as can be expected from dimensional arguments, this param¬ 
eter is a function of the collapse redshift. We compare the cir¬ 
cular velocity of the final NFW and corresponding a model 
profiles. The values are similar at low circular velocity (low 
c or a) and diverge for high circular velocity, but are close to 
one another for the relevant range of circular velocity values 
seen in dark matter haloes. Thus v‘^ax/{GM 2 oo/R 20 o) could 
also potentially be used as a parameter in a halo matching 
scheme, and would have the benefit of not being highly de¬ 
pendent on the fitting function chosen for the dark matter 
halo profile. 



0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 

X’max/{GM200/R200) 


Figure 3. The normalized squared maximum circular velocity, 
/{GM200/R200 ) I versus the normalized excess central po¬ 
tential, / (G M 20 Q/ R200) t at z = 0 for an NFW profile with 

varying c and an a model with varying inner slope a. Unlike the 
maximum circular velocity shown in Figure (21 the central poten¬ 
tial is significantly different for the two models. 


In Figure O we show the v^ax/ {GM200/R200) versus a 
related parameter, the excess normalized central potential 
5 ^/{GM200/R200), for both NFW profiles and a models. 
While the values of v^a,x / (GM200 / R200 ) are similar for both 
models at fixed mass and energy, the values of the central po¬ 
tential are signihcantly different for the two profiles shapes, 
implying that the maximum circular velocity is a superior 
“common” parameter between the two types of models to 
use to predict the initial halo collapse redshift. 

For reference, we show contours of constant '^max in 
Figure [4l in the top panel for varying M200 and c for NFW 
profiles, and in the bottom panel for varying M200 and a for 
a models, over the range of physical interest. 


4 COMPARISON TO PREVIOUS WORK 


Previous papers have investigated parametrizations of 
the time at which a halo formed, usually based on 
the r esults of dark matt er si mulations. Such pape rs in¬ 
clude |Bunoc^|eF^i]jl22Qll) and IWechsler et al.l (I2OO2I I. The 


model of Wec hsler et ahf 1120021 improved upon that of 
[Bullock et al.l 1I2OOI ) , but obtained similar results. 

IWechsler et al.l 1I2OO2II examined the mass accretion his¬ 
tories of the most massive progenitors of individual haloes 
from dark matter simulations for z < 7, dehning the hal o 
mass to be the mass within Avir llBrvan fc Normarilll998ll . 
IWechsler et al.l (I2OO2I ') found that the mass accretion histo¬ 
ries of haloes at each timestep, despite halo mergers in the 
simulation, can generally be fit well by a simple analytic 
form: 


M{a) = Mo exp |^—acS" “ l)] ■ 


( 20 ) 


Here a is the scale factor, ao is a reference scale factor at 
which the mass of the halo is Mo, S is an arbitrarily chosen 
constant, and is taken to be the “formation scale factor” 
of the halo, given the choice of S. This form is self-consistent 
for different choices of oq. Taking the formation scale factor 
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Table 1. Central collapse redshifts for an initial Gaussian profile at Zi = 1000. This table shows the central collapse redshifts of 
initial Gaussians matched to NFW profiles with varying M 2 qq and c . Bold numbers indicate the one-sigma scatter around the z = 0 
mass-concentration relation for NFWs from lPiemer fc Kravtso'^ ll2014li . 



l.OOe + lO 

2.51e+10 

6.31e+10 

1.58e+ll 

^200 [-^©1 

3.98e+ll 

l.OOe+12 

2.51e+12 

6.31e+12 

1.58e+13 

3.0 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

1.56e+00 

4.0 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e+00 

2.11e-h00 

5.0 

2.73e+00 

2.73e+00 

2.73e+00 

2.73e+00 

2.73e+00 

2.73e+00 

2.73eS-00 

2.73e-h00 

2.73e-h00 

6.0 

3.40e+00 

3.40e+00 

3.40e+00 

3.40e+00 

3.40e-h00 

3.40e-h00 

3.40e-{-00 

3.40e-h00 

3.40e-h00 

7.0 

4.15e+00 

4.15e+00 

4.15e+00 

4.15e-hOO 

4.15e-f-00 

4.15e-h00 

4.15ef-00 

4.15e-h00 

4.1Se-hOO 

8.0 

5.00e+00 

5.00e+00 

5.00e-h00 

5.00e-f-00 

S.OOe-hOO 

S.OOe-hOO 

S.OOe-hOO 

S.OOe-hOO 

5.00e+00 

9.0 

5.94e+00 

5.94e-fOO 

5.94e-hOO 

5.94e-f00 

5.94e-h00 

5.94e + 00 

5.94e-h00 

5.94e+00 

5.94e+00 

10.0 

7.04e-l-00 

7.04e-f00 

7.04e-l-00 

7.04e-h00 

7.04e-h00 

7.04e + 00 

7.04e+00 

7.04e+00 

7.04e+00 

11.0 

8.29e-h00 

8.29e-h00 

8.29e-h00 

8.29e-h00 

8.29e-h00 

8.29e-h00 

8.29e+00 

8.29e+00 

8.29e+00 

12.0 

9.76e-f00 

9.76e-h00 

9.76e-f00 

9.76e-f00 

9.76e-h00 

9.76e+00 

9.76e+00 

9.76e+00 

9.76e+00 

13.0 

l.lSe-j-Ol 

l.lSe-fOl 

1.15e-f01 

l.lSePOl 

1.15e+01 

1.15e+01 

1.15e+01 

1.15e+01 

1.15e+01 

14.0 

1.36e-j-01 

1.36e-h01 

1.36e-h01 

1.36e+01 

1.36e+01 

1.36e+01 

1.36e+01 

1.36e+01 

1.36e+01 

15.0 

1.63e-h01 

1.63e-h01 

1.63e-j-01 

1.63e+01 

1.63e+01 

1.63e+01 

1.63e+01 

1.63e+01 

1.63e+01 

16.0 

1.97e-h01 

1.97e-h01 

1.976+01 

1.97e+01 

1.97e+01 

1.97e+01 

1.97e+01 

1.97e+01 

1.97e+01 

17.0 

2.42e-f-01 

2 . 42 e-h 01 

2.42e+01 

2.42e+01 

2.42e+01 

2.42e+01 

2.42e+01 

2.42e+01 

2.42e+01 

18.0 

3.03e-f01 

3.03e+01 

3.03e+01 

3.03e+01 

3.03e+01 

3.03e+01 

3.03e+01 

3.03e+01 

3.03e+01 

19.0 

3.93e-j-01 

3.93e+01 

3.93e+01 

3.93e+01 

3.93e+01 

3.93e+01 

3.93e+01 

3.93e+01 

3.93e+01 

20.0 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 

5.31e+01 


Table 2. Half-mass collapse redshifts for an initial Gaussian profile at Zi = 1000. Same as Tabled but for the half-mass collapse redshift. 



l.OOe+lO 

2.51e+10 

6.31e+10 

1.58e+ll 

^200 [-^©1 

3.98e+ll 

l.OOe+12 

2.51e+12 

6.31e+12 

1.58e+13 

3.0 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.87e-01 

4.0 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

6.07e-01 

5.0 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

7.18e-01 

6.0 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

8.19e-01 

7.0 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

9.12e-01 

8.0 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.99e-01 

9.0 

1.08e+00 

1.08e-h00 

1.08e-h00 

1.08e-h00 

1.08e-h00 

1.08e-h00 

1.08e-h00 

1.08e+00 

1.08e+00 

10.0 

l.lSe-hOO 

l.lSe-hOO 

l.lSe-hOO 

l.lSe-hOO 

l.lSe-hOO 

l.lSe-hOO 

1.15e+00 

1.15e+00 

1.15e+00 

11.0 

1.22e-h00 

1.22e-h00 

1.22e-h00 

1.22e-h00 

1.22e-h00 

1.22e-h00 

1.22e+00 

1.22e+00 

1.22e+00 

12.0 

1.28e-h00 

1.28e-h00 

1.28e-h00 

1.28e-h00 

1.28e-h00 

1.28e+00 

1.28e+00 

1.28e+00 

1.28e+00 

13.0 

1.34e-hOO 

1.34e-hOO 

1.34e-hOO 

1.34e + 00 

1.34e+00 

1.34e+00 

1.34e+00 

1.34e+00 

1.34e+00 

14.0 

1.38e-h00 

1.38e-h00 

1.38e-h00 

1.38e+00 

1.38e+00 

1.38e+00 

1.38e+00 

1.38e+00 

1.38e+00 

15.0 

1 . 42 e-hOO 

1 . 42 e-h 00 

1 . 42 e-h 00 

1.42e+00 

1.42e+00 

1.42e+00 

1.42e+00 

1.42e+00 

1.42e+00 

16.0 

1.46e-hOO 

1.46e-hOO 

1.46e+00 

1.46e+00 

1.46e+00 

1.46e+00 

1.46e+00 

1.46e+00 

1.46e+00 

17.0 

1.48e-hOO 

1.48e-h00 

1.48e+00 

1.48e+00 

1.48e+00 

1.48e+00 

1.48e+00 

1.48e+00 

1.48e+00 

18.0 

1.49e-hOO 

1.49e+00 

1.49e+00 

1.49e+00 

1.49e+00 

1.49e+00 

1.49e+00 

1.49e+00 

1.49e+00 

19.0 

l.SOe-hOO 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

20.0 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 

1.50e+00 


ttc of the halo as defined above with S chosen to be 2, the 
authors assign a formation time to each halo and find that 
the concentration of each halo at some reference scale factor 
ao is related to the formation time via 


Cvir -— 4.1(2 o/o 


( 21 ) 


where the constant 4.1 is the concentration of haloes “form¬ 
ing” at the present day, given the choice of S' = 2. 

A more recent paper lo oking at ma.ss accr etion histories 
in dark matter simulations is IZhao et al.l (|2009l l. The authors 
claim to obtain more accurate fit s to halo accretion hist ories 
using a more c o mplex model than Wechsler et al.l (l2002f l and 
iBullock et al.l (1200 il l. Unlike IWechsler et al.l ^ 2(302h . they 
find that mass accretion histo ries follow a powe r-law form 
as opposed to an exponential. IZhao et al.l (l2009tl find that 
the concentration of haloes at some observation time t is a 
function of the time at which their mass reached 4% of the 
mass at the observation time, given by 


C=[4« + (t/to.04)®-1'/® 


( 22 ) 


where to,o 4 is the time at which 4% of the mass at t was 


reached. Note that halo concentrations in this model cannot 
be less than 4. 

Due to the fact that we begin with arbitrary initial 
conditions and not cosmological ones, it is difficult to com- 
pare our results with those of IWechsler et al.l ll2002ll and 
IZhao et ^ ll2009tl . However, both our model and the results 
obtained by these two papers find a monotonic relation be¬ 
tween the collapse time or “formation time” of haloes and 
their concentrations at fixed mass at any given time, where a 
higher concentration implies an earlier formation time. This 
is crucial as abundance matching-type methods require a 
monotonic relation between the parameters being matched. 


5 FUTURE WORK 

Using a combination of two halo parameters, one can per¬ 
form two-parameter abundance matching to two observ¬ 
able parameters of galaxies. As shown above, parameters 
that correlate well with the physically relevant parame¬ 
ter of halo collapse time include c or a, or alternatively 
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12.0 12.5 13.0 13.5 14.0 14.5 15.0 


log (A/200 [Ai'o]) 

Figure 4. Contours of constant Vmax. for varying M 200 and c for 
NFW profiles in the top panel, and the same for varying M 200 
and inner slope for a models in the bottom panel. 


Pmax!{GM 200 /R 20 o)- These could be combined with the pa¬ 
rameter M 200 used in standard abundance matching, al¬ 
though these are not the only possible choices for matching. 

Because we only consider two parameters, and standard 
abundance matching is most effective for massive galaocies, 
we plan to focus on two-parameter abundance matching to 
the Fundamental Plane of massive early type galaxies in 
future work. The Fundamental Plane is an observed rela¬ 
tionship between three parameters: the effective radius Re, 
the velocity dispersion cr, and the surface brightness within 
Re, le (or alternatively, the luminosity within Re). Elliptical 
galaxies occupy a plane in the space of these three parame¬ 
ters: 

log /?e = a log cr -I- fe log le + Cz, (23) 

where a and b are constants and Cz is a redshift-dependent 
zero-point. Assuming that the mass-to-light ratio of elliptical 
galaxies is roughly independent of stellar mass, and that the 
galaxies are fully virialized, one would expect from the virial 
theorem that the fundamental plane would follow Re oc 
Ple^. However, the observed fu ndamental plane has a tilt 
with respect to this relation (e.g.. lBezanson et al 11201311 . Al¬ 
though the FP has some thickness and i s not completely two- 
dimensional, this thic kness is small llGraves et akl l2009bl : 
[Graves &: Fabeil l20ld l . W hile we plan to focus on early- 
type galaxies, interestingly, iBezanson et al.l ll2014ll find that 


early- and late-type galaxies fall on the same mass funda¬ 
mental plane. 

Additionally, other galactic parameters have been found 
to be strongly correlated with parameters of the Funda¬ 
mental Plane. These include galaxy color, stellar mass-to- 
light ratio, and mean st ellar population age a n d metallicity 
( [Graves et al.l I2009al f9: [Graves fc Fabeil 1(^101 : IPorter et al.l 
[2014h . Thus matching to the FP would take into account 
the variance in these parameters as well, potentially allow¬ 
ing a near-complete prediction of the properties of a massive 
galaxy residing in a halo based off the halo’s properties, and 
creating a better connection between the properties of the 
halo (including formation epoch) and the properties of the 
galaxy. 
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APPENDIX 

In tl2.1l we show the equations determining the evolution of 
mass shells of an overdensity in a ACDM universe, where 
each shell has initial condition r = 0 at t = 0. In a ACDM 
universe, mass shells with overdensity S below some value 
will never turn around. For the case presented in am there 
is no simple analytic formnla for this overdensity for all times 
(although it can be found numerically from the equations 
given therein), but there is one as the initial time approaches 
0, given by Eqn. 1131 However, for the case in which all mass 
shells are taken to be initially on the Hubble flow (i.e., ri = 
Hi[ti)ri for a chosen initial time ti), there is an analytic 
solution for the minimnm overdensity to turn around for all 


initial times. Here we present a short derivation of this value 
and compare to the solution for initial conditions r = 0 at 
t = 0 at early times. 

Once again, we begin with the equation of motion for a 
shell containing mass M: 


Pr 


GM 2 

—^—h Da,oJ7o*'- 


r 

Integrated once, we obtain 

1 .2 1 .2 _ GM GM IlA.oiJo 

2^ ~ 2 ^'“ ~ ~~r ~ 2 

Here our initial condition is 


2 Q,A,otiQ 2 
r- - - Ta . 


(24) 


(25) 


.2 2 2G{M) rrS 2 

Vi = HiVi =-h ilA,oHori, 


(M) = -nGp{ti)rf. 

(26) 


Combined with Eqn. 

2GM 


we obtain 


.2 

r = 


+ DA,o77oV-Mill^M). 


(27) 


Now we define 


/ - d _ d 


(28) 


giving 

dx 

dr 


8nGp M nA,oH§ 


{M)x 


H? 


2 SttGp M — (M) 
X 


3Hf (M) 


Noting that 


l-bJ = 


M 


StvGp 


— 


(M) ’ 3Hf 
and that for a flat universe 
2 


H? 


(29) 


= Da.i (30) 


m,^ + Da,; = 1, we obtain 

^^d(l J” d') , ^ 2 I", 

- -}-ilA.iX — (31) 


For any shell that turns around, (dx/dr)^ = 0 for some 
maximal value of x- Therefore we want to find the smallest 
5 such that there exists a positive x where 


X 


c- , 

-dx + 


-(1 + S ) = 0 . 


Da,! ' Da,* ^ ^ ^ ^ 

We also require that this cubic to be positive for x between 
0 and the maximum value of x in order for dx/dr to be 
a real number. This cubic has inflection points at x± = 
±(Dm,*d/3DA,*)^^^; for the previously mentioned conditions 
to hold, it requires that the value of the cubic at x+ be 
negative or 0 : 

3/2 


/ ^m,i 


f r ^ 

k DA.i 3^ 

1 31/2 1 

^DA,i J 




Da,* 


This reduces to simply 


> 


27 Dy. 


(1-f 5) < 0. (33) 


(34) 


(l-bd )2 - 4 Da,! ' 

As stated, this is the exact solution for an initial profile 
that is assumed to be on the Hubble flow at time ti. This 
can be compared to Equation 1131 which is the solution for 
5 with initial conditions r = 0 at t = 0 for all shells as 
Zi oo. Here the boundary between collapsing and forever 
expanding shells is given by (" = 0.5, by the definition of 
This gives for the minimal overdensity at early times 

1/3 

; (35) 


J = -^2^/3 
10 
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whereas at early times for the solution on the Hubble flow 
we find 




2 ^ 



(36) 


It can be seen that the derived S differs by a value of 
3/5. This is the difference in overdensity expected for a 
matter-dominated universe between considering a hypersur- 
face of constant Hj and one of constant ti, as derived in 
I Gunn fc Gotti lll972l l. This ratio between the results will not 
hold at small Zi due to the effect of A, but for early times 
the universe is increasingly matter dominated and so the 
two solutions are the same modulo this factor of 3/5. 
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